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Abstract. A potential model description of heavy quarkonium can be rigorously deduced from 
QCD under some circumstances. The potentials can be unambiguously related with Wilson loops 
with gluonic insertions, the spectral decomposition of which is a function of the spectrum and 
matrix elements solution of the static limit of NRQCD. This spectrum is nothing but the static 
singlet potential and the hybrid potentials (which correspond to the gluonic excitations). We will 
quantitatively show that the latter unambiguously relate to the gluelumps at short distances using 
effective field theories. 
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INTRODUCTION 

Potentials appear in a natural way in the description of systems with slow moving parti- 
cles like in heavy quarkonium. It is now possible 01 12] to obtain a potential model-like 
description of the heavy quarkonium dynamics from first principles using effective field 
theories under some circumstances (for a review see (3D). The study of these potentials 
(the static [4] and the relativistic |5i la LZD corrections) becomes then crucial as they dic- 
tate the dynamic of the heavy quarkoniurrQ. They can be obtained from the solution 
(spectrum and matrix elements) of NRQCD in the static limit. The spectrum of NRQCD 
directly provides the static singlet potential and the hybrid potentials. Therefore it is 
interesting to study them on their own. In particular, the study of the hybrid potentials 
may become important in order to discern whether the recently found resonances near 
threshold can be understood as hybrid states (see [8] for a review on the different possi- 
ble interpretation of these new states). 

In this review we closely follow [9]. We will focus on the short distance limit of the 
hybrid potentials, specially on the aspects that can be fixed on symmetry arguments only, 
and we will quantitatively relate them with the physics of gluelumps. 



One may wonder why bother to compute these potentials in the lattice, as one may try to do a direct 
computation of the heavy quarkonium properties in the lattice. Leaving aside the technical difficulties of 
these direct computations, in the best of the worlds they only provide with a very limited information on 
the properties of the heavy quarkonium: mainly a few masses (for the ground state and low excitations) 
and, maybe, some inclusive decays. On the other hand, with the potentials it is possible to obtain a detailed 
information of the shape of the heavy quarkonium, one can then compute the complete spectrum, all the 
inclusive decays and also opens the possibility to consider differential decay rates, which are far beyond 
the possibility of present direct lattice simulations. 



STATIC LIMIT OF PNRQCD 



The static limit of NRQCD at short distances can be studied with the static version of 
pNRQCD. In this limit new symmetries arise. 

The p NRQ CD Lagrangian at leading order in 1 jm and in the multipole expansion 
reads IKALLil], 

^pNRQCD = Jd\d i RTx S f (id -V s )S+O* (iD -V o )O - J d 3 R^ v F^ Va + 0(r). 

(1) 

All the gauge fields in Eq. (OQ) are evaluated in R and t, in particular F^ Va = F^ lVa (R 1 t) 
and iDqO = idoO — g[Ao(R,?),0], The singlet and octet potentials Vu i = s,o are to be 
regarded as matching coefficients, which depend on the scale v us separating soft gluons 
from ultrasoft ones. In the static limit "soft" energies are of 0(1 /r) and "ultrasoft" 
energies are of 0(a s /r). Notice that the hard scale, m, plays no role in this limit. The 
only assumption made so far concerns the size of r, i.e. l/r> Aqcd> such that the 
potentials can be computed in perturbation theory. Also note that throughout this paper 
we will adopt a Minkowski space-time notation. 

The spectrum of the singlet state reads, 

E s (r)=2m os + V s (r)+0(r 2 ), (2) 

where mos denotes an on-shell (OS) mass. One would normally apply pNRQCD to 
quarkonia and in this case mos represents the heavy quark pole mass. For the static 
hybrids, the spectrum reads 

E H {r) = 2m os + V {r) + A° s + 0(r 2 ) , (3) 

where 

A° s = r lhn^ln(// fl (r/2)<^(r/2,-r/2)// fc (-r/2)). (4) 

0(772,-772) = 0(772, R, -r/2, R) 

= ?Qxp^-igj^dtA (R,t)j, (5) 

denotes the Schwinger line in the adjoint representation and H represents some gluonic 
field, for examples see Table [TJ 

Eq. © allows us to relate the energies of the static hybrids Eh to the energies of the 
gluelumps, 

A° s = [E H (r)-E s (r)]-[V o (r)-V s (r)] + 0(r 2 ). (6) 

This equation encapsulates one of the central ideas of this paper. The combination 
Eh — E s is renormalon-free in perturbation theory [up to possible 0(r 2 ) effects], and 
can be calculated unambiguously non-perturbatively: the ultraviolet (UV) renormalons 
related to the infrared (IR) renormalons of twice the pole mass cancel each other. 
However, Ah contains an UV renormalon that corresponds to the leading IR renormalon 
ofV . 



Symmetries of hybrid potentials and gluelumps 



The spectrum of open QCD string states can be completely classified by the quantum 
numbers associated with the underlying symmetry group, up to radial excitations. In 
this case, these are the distance between the endpoints, the gauge group representation 
under which these endpoints transform (in what follows we consider the fundamental 
representation), and the symmetry group of cylindrical rotations with reflections D^. 
The irreducible representations of the latter group are conventionally labelled by the 
spin along the axis, A, where E,n, A refer to A = 0, 1, 2, respectively, with a subscript 
7] = g for gerade (even) PC = + or r\ = u for ungerade (odd) PC = — transforma- 
tion properties. All A > 1 representations are two-dimensional. The one-dimensional E 
representations have, in addition to the 77 quantum number, a o v parity with respect to 
reflections on a plane that includes the two endpoints. This is reflected in an additional 
± superscript. The state associated with the static singlet potential transforms according 
to the representation E+ while the two lowest lying hybrid potentials are within the U u 
and E~ representations, respectively. 

In contrast, point-like QCD states are characterised by the J of the usual 0(3) ® ^ 
rotation group as well as by the gauge group representation of the source. In the pure 
gauge sector, gauge invariance requires this representation to have vanishing triality, 
such that the source can be screened to a singlet by the glue. States created by operators 
in the singlet representation are known as glueballs, octet states as gluelumps. In contrast 
to gluelump states, where the octet source propagates through the gluonic background, 
the normalization of glueball states with respect to the vacuum energy is unambiguous. 

Since C 0(3) <g> < €, in the limit r — > certain hybrid levels must become degen- 
erate. For instance, in this limit, the E~ state corresponds to a J = l + ~ state with 
J z = while the II [( doublet corresponds to its J z = ±1 partners. The gauge transfor- 
mation property of the hybrid potential creation operator will also change in this limit, 
3 ® 3* = 1 © 8, such that hybrids will either approach gluelumps [cf. Eq. ©] or glue- 
balls, in an appropriate normalization. In the case of glueballs the correct normalization 
can be obtained by considering the difference En{r) —E s (r) from which the pole mass 
cancels. We will discuss the situation with respect to gluelumps in detail below. 

In perturbation theory, the ground state potential corresponds to the singlet potential 
while hybrid potentials will have the perturbative expansion of the octet potential, which 



We would like to establish if lattice data on hybrid potentials reproduces the degen- 
eracies expected from the above discussion in the short distance region. In the limit 
r — > 0, any given A > 1 hybrid potential can be subduced from any J PC state with 
J > A and PC = + for 77 = g or PC = — for 77 = u representations. For instance the 
U u is embedded in 1 + ~,1 h ,2 +_ ,2 h ,---. The situation is somewhat different for 
A = states, which have the additional o v parity: the E+ representation can be ob- 
tained from 0++ y l—,2 ++ ,--; E„ from 0" ,1 ++ ,---, E+ from + ~, !-+,■•• and E" 




Hybrid and gluelump mass splittings 



TABLE 1. Expected degeneracies of hybrid potentials at short distance, based 
on the level ordering of the gluelump spectrum. Note that if the 3 H gluelump 
turned out to be lighter than the 2 + ~ then the £~ ,Yl' u , A„,4>„ potentials would 
approach the 3 H state while the E+,njJ,Aj, potentials would approach the 2 H 
instead. 
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from h , 1 H , • • •. We list all combinations of interest to us in Table [Q The ordering of 
low lying gluelumps has been established in Ref. dl2D and reads with increasing mass: 
l+~, 1 ,2 ,2+~,3+~,0 ++ ,4 , with a 3 state in the region of the 4 and 
1~ + . The 2 H and 3 +_ as well as the 4 and 1~ + states are degenerate within present 
statistical uncertainties. The continuum limit gluelump masses are displayed as circles 
at the left of Fig. [T] where we have added the (arbitrary) overall constant 2.26/ro to the 
gluelump splittings to match the Jiybrid potentials. 

Juge, Kuti and Morningstar lfl3[ [Hi have comprehensively determined the spectrum 
of hybrid potentials. We convert their data, computed at their smallest lattice spacing 
a a ~ 0.2 fm, into units of tq w 0.5 fm. Since the results have been obtained with an 
improved action and on anisotropic lattices with I A, one might expect lattice 

artifacts to be small, at least for the lower lying potentials. Hence we compare these data, 
normalized to E £ + (ro), with the continuum expectations of the gluelumps [12]. The full 
lines are cubic splines to guide the eye while the dashed lines indicate the gluelumps 
towards which we would expect the respective potentials to converge. 

The first 7 hybrid potentials are compatible with the degeneracies suggested by Ta- 
bled] The next state is trickier since it is not clear whether 2 H or 3 H is lighter. In the 
figure we depict the case for a light 2^ . This would mean that (E+,n^,A u ) approach 
the 2 + ~ while (L~', IT", A' u , 4> M ) approach the 3 + ~. Note that of the latter four potentials 
only data for II" and 4> ;( are available. Also note that the continuum states IT^, IT" and 4> ;( 
are all obtained from the same E u lattice representation. For the purpose of the figure we 
make an arbitrary choice to distribute the former three states among the E' U ,E„ and E'// 
lattice potentials. To firmly establish their ordering one would have to investigate radial 
excitations in additional lattice hybrid channels and/or clarify the gluelump spectrum in 
more detail. Should the 2 + ~ and 3^ hybrid levels be inverted then ,Tl' u ,A u ,Q> u ) 
will converge to the 3^ while (E+,n",A^) will approach the 2^ . We note that the 



The splittings between all states with respect to the l" 1 ground state have been extrapolated to the 
continuum limit in Ref. uA\ and we add our own extrapolations for the 4 and 1 h states to these, 
based on the tables of this reference. 
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FIGURE 1. Different hybrid potentials lfl4ll at a lattice spacing a a w 0.2 fm ss 0.4 rp where ro ~ 0.5 fm, 
in comparison with the gluelump spectrum, extrapolated to the continuum limit 111211 (circles, left-most 
data points). The gluelump spectrum has been shifted by an arbitrary constant to adjust the l + ~ state 
with the n„ and E,7 potentials at short distance. In addition, we include the sum of the ground state (E+) 
potential and the scalar glueball mass m Q ++. The lines are drawn to guide the eye. 

ordering of the hybrid potentials, with a low E+, makes the first interpretation more 
suggestive. 

Finally the E+" potential seems to head towards the ++ gluelump but suddenly turns 
downward, approaching the (lighter) sum of ground state potential and scalar glueball 
instead. The latter type of decay will eventually happen for all lattice potentials but only 
at extremely short distances. We also remark that all potentials will diverge as r — ► 0. 
This does not affect our comparison with the gluelump results, since we have normalized 
them to the II M /E~ potentials at the shortest distance available. (The gluelump values 
are plotted at r = to simplify the figure.) 

On a qualitative level the short-distance data are very consistent with the expected 
degeneracies. From the figure we see that at r ~ 2tq lfm the spectrum of hybrid 
potentials displays the equi-distant band structure one would qualitatively expect from a 
string picture. Clearly this region, as well as the cross-over region to the short-distance 
behaviour tq < r <2tq, cannot be expected to be within the perturbative domain: at best 
one can possibly imagine perturbation theory to be valid for the left-most 2 data points. 



With the exception of the Tl u , Tl' u and 4>„ potentials there are also no clear signs for 
the onset of the short distance 1/r behaviour with a positive coefficient as expected from 
perturbation theory. Furthermore, most of the gaps within multiplets of hybrid potentials, 
that are to leading order indicative of the size of the non-perturbative r 2 term, are still 
quite significant, even at r = 0.4 ro; for instance the difference between the £~ and U u 
potentials at this smallest distance is about 0.28 r Q 1 « 1 10 MeV. 



The difference between the n M and H u hybrids 

From the above considerations it is clear that for a more quantitative study we need 
lattice data at shorter distances, which have been obtained for the lowest two gluonic 
excitations, U u and E~ in Ref. Q9[| . We display their differences in the continuum limit 
in Fig. [2l We see how these approach zero at small r, as expected from the short distance 
expansion. pNRQCD predicts that the next effects should be of 0(r 2 ) (and renormalon- 
free). In fact, we can fit the lattice data rather well with a Ejj u — E L = A u _ £ -r 2 ansatz 
for short distances, with slope (see Fig. [2]), 

^n„-z l7 =0.92lS;gr 3 , (7) 

where the error is purely statistical (lattice). This fit has been done using points 
r < 0.5 ro. By increasing the fit range to r < 0.8ro the following result is obtained, 

A u u -^ = (0-83 ± 0.29) r - 3 , (8) 

indicating stability of the result Eq. © above. 

In order to estimate systematic errors one can add a quartic term: br 4 (only even 
powers of r appear in the multipole expansion of this quantity). If the result is stable, our 
determination of A n should not change much. Actually this is what happens. If we 

fit up to r < 0.5 ro, we obtain the central value An„-£,/o = 0.93 with a very small quartic 

coefficient, br^ = —0.05. If we increase the range to r < 0.8 ro, we obtain the same 

central value, A n = 0.93, but with a slightly bigger quartic term, br^ = —0.18. 

Introducing the quartic term enhances the stability of A n £ under variations of the 
fit range. From this discussion we conclude that the systematic error is negligible, in 
comparison to the error displayed in our result Eq. ©. 

We remark that within the framework of static pNRQCD and to second order in the 
multipole expansion, one can relate the slope A n _ L - to gluonic correlators of QCD. 



CONCLUSIONS 

In conclusion, we have shown that the short distance behavior of the hybrid potentials 
is consistent with perturbation theory, the operator product expansion (effective field 
theories), and with the description of the leading non-perturbative effects in this regime 
in terms of the gluelump masses. For more details see 
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FIGURE 2. Splitting between the E,7 and the n„ potentials, extrapolated to the continuum limit, and the 
comparison with a quadratic fit to the r < 0.5 ro data points (r^ 1 0.4 GeV). The big circles correspond 
to the data of Juge et al. 11411 . obtained at finite lattice spacing a a w 0.39 ro. The errors in this case are 
smaller than the symbols. 
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